LA-UR-09-03593 
NUHEP-TH/09-08 



Pseudo-Dirac Neutrinos in the New Standard Model 



Andre de Gouvea and Wei-Chih Huang 
Northwestern University, Department of Physics & Astronomy, 2145 Sheridan Road, Evanston, IL 60208, USA 

James Jenkins 
Elementary Particles and Field Theory Group, MS B285, 
Los Alamos National Laboratory, Los Alamos, NM 87545, USA and 
Northwestern University, Department of Physics & Astronomy, 2145 Sheridan Road, Evanston, IL 60208, USA 

The addition of gauge singlet fermions to the Standard Model Lagrangian renders the neutrinos 
massive and allows one to explain all that is experimentally known about neutrino masses and 
lepton mixing by varying the values of the Majorana mass parameters M for the gauge singlets and 
the neutrino Yukawa couplings A. Here we explore the region of parameter space where M values 
are much smaller than the neutrino Dirac masses \v. In this region, neutrinos are pseudo-Dirac 
fermions. We find that current solar data constrain M values to be less than at least 10~^ eV, and 
discuss the sensitivity of future experiments to tiny gauge singlet fermion masses. We also discuss a 
useful basis for analyzing pseudo-Dirac neutrino mixing effects. In particular, we identify a simple 
relationship between elements of M and the induced enlarged mixing matrix and new mass-squared 
differences. These allow one to directly relate bounds on the new mass-squared difi'erences to bounds 
on the singlet fermion Majorana masses. 



Nonzero neutrino masses reveal that the minimum standard model needs to be modified. Current data (mostly 
from solar, atmospheric, reactor and accelerator neutrino oscillation experiments), however, provide only minimal 
insight as to how the standard model ought to be extended. Several completely different new physics scenarios can 
be constructed and all safely agree with observations, i.e., all lead to small neutrino masses and non-trivial mixing 
among the three so-called active neutrino weak eigcnstates, Vf., v^j,, and Vr- Additions to the standard model that fit 
the data include right-handed neutrinos, SU (2) l triplet fermions or scalars, lepto-quarks, etc. Recent summaries of 
the current experimental and theoretical situation can be found, for example, in Q, S i, i, S, H H ■ 

One version for the new standard model, :/SM, consists of the minimum standard model augmented by a few (at 
least two) SU{2)l x U{1)y gauge singlet Weyl fermions N, normally referred to as right-handed neutrinos. At the 
renormalizablc level, the most general I'SM Lagrangian consistent with SU{3)c x SU{2)l x U{1)y gauge invariance is 



where £oid is the minimal standard model Lagrangian, L are the SU{2)l lepton doublet fields L = (y t) and H is the 
standard model Higgs doublet field. Aq^ are neutrino Yukawa couplings and Afy are Majorana masses for the N fields. 
Note that M is a symmetric matrix: = Mji. a = e, ^, r, i = 1, . . . n, and n > 2 is the number of right-handed 
neutrino fields. After clcctrowcak symmetry breaking, Eq. p.l[) describes 3 -I- n Majorana neutral fermions, referred 
to as neutrinos. In general, all neutrino mass eigenstates vi,V2t ■ ■ ^ (with mass mi, m2, . . . , TO3+„) are linear 
combinations of the active neutrinos plus the n right-handed neutrinos. Since these singlet fermions do not couple to 
any of the standard model gauge bosons, they will also be referred to as sterile neutrinos. 

The phenomenology associated to Eq. (jl.l[) depends dramatically on the values of the i^SM parameters A and M 
(see, for example, Q). Even after the current data are taken into account, the spectrum of possibilities remains vast. 
When the matrix M vanishes exactly, the six* neutral states fuse into three Dirac fermions with masses proportional 
to the (square-roots of the) eigenvalues of A^A. In this case, the classical global symmetry of £^sm is enhanced and 
includes the lepton-numbcr symmetry U{\)l (and its non-anomalous cousin U{1)b-l)- This indicates that AI can be 
interpreted as a symmetry breaking parameter (when A 7^ 0) so that (i) quantum corrections to M are proportional 
to M itself and hence (ii) any value of AI is technically natural. 

When M 7^ 0, one can identify three qualitatively distinct regions of the A, AI parameter space. In the region where 
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* For the remainder of this section, we restrict our discussion to n = 3 right-handed neutrinos. 
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M ^ Aw,''' where v is the vacuum expectation value of the neutral component of H, the six neutrino masses "split" 
into three lighter, mostly active states with masses generically of order X^v'^lM and three heavier, mostly sterile states 
with masses M. This phenomenon is referred to as the seesaw mechanism and it was argued in Q that all values 
of Af > 1 eV are phenomenologically allowed. For small enough AI (in general M much smaller than the weak scale), 
the z^SM can be tested experimentally since the mostly sterile, heavier states are accessible at different facilities (see, 
for exarnple, [HI, ) • It has also been pointed out that the mostly heavy states may qualify as viable warm dark 
matter [l^. Finally, for extraordinary choices of A (for recent discussions see [13, Ell): rnostly sterile neutrinos with 
weak-scale masses can be detected at high energy collider experiments. 

The other two distinct regions of the A, M parameter space are characterized by M ^ \v and M <C Au. In 
the former, all six neutrino masses are of the same order, and all six neutrino mass eigenstates are characterized 
by "homogeneous" mixtures of active and sterile flavors. Such a scenario is both hard to study quantitatively and 
severely constrained by current solar and atmospheric neutrino data. In this paper, we explore the latter possibility: 
right-handed neutrino Majorana masses M much smaller than the so-called Dirac neutrino masses Xv. Under these 
circumstances, neutrinos are pseudo-Dirac fermions [3, [13, [a. 

Pseudo-Dirac neutrinos are Majorana neutrinos made up of roughly fifty-fifty mixtures of active and sterile neutri- 
nos, and come in quasi-degenerate pairs. These will be properly defined in Sec. [Til In the limit M 0, neutrinos are 
Dirac fermions and hence can accommodate all current experimental data (as discussed above) so our main goal is 
to estimate an upper bound for AI . Given our current understanding of neutrinos, the most stringent constraints on 
very small M values are provided by solar data. The data, the estimation procedure, and our results are presented in 
Sec. mil Other constraints and predictions, as well as expectations for the future, are discussed in Sec. |TVl A sum- 
mary of what is currently known about pseudo-Dirac z/SM neutrinos and the i/SM Lagrangian (Eq. (|I.1[) ) in general 
is presented in Sec. |Vl 



II. "ANTI-SEESAW:" MASSES AND MIXING 



In the case of 7i = 3 right-handed neutrinos, after electroweak symmetry breaking, the 6x6 Majorana neutrino 
mass matrix is 

- ( % T ) , (n.i) 



where O3 stands for the 3x3 zero matrix, 



V^e^Vn, (II.2) 



is the symmetric 3x3 Majorana right-handed neutrino mass matrix, and wc choose the weak basis where the 3x3 
Dirac neutrino mass matrix is written as 

m = U*m^. (II.3) 

This can be achieved by redefining the right-handed neutrino fields N. Above, the superscript D indicates a diagonal 
matrix while and U are unitary matrices. 

Since we are interested In the limit ^ - Dirac masses much larger than Majorana right-handed neutrino 
masses - we can write 

where I3 is the 3x3 unit matrix, and the elements of S are small (this will be shown a posteriori) so that 



I3 S \ f l3 -S 

-jt I3 i Ut 13 



(II.5) 



t Beware of the abuse of notation and keep in mind that A and M are matrices. For a more detailed discussion of this so-called seesaw 
limit, see, for example, Q. 
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is a diagonal dimcnsionlcss matrix of small numbers. Eq. (jll.4[) is satisfied, at leading order, if 



(II.6) 



where e is a dimensionless matrix of small numbers whose diagonal elements vanish. Both and e are functions of 
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and ~ —em^ . 



(11.7) 
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Note that and e do not commute and that Cm is indeed a symmetric matrix with units of mass. Throughout 
we will ignore the possibility, ruled out by data, that different diagonal entries of are identical. According to 
Eq. (|II.7[) . in the weak basis of choice and at leading order, the diagonal elements of e„i, proportional to e^, determine 
the mass-squared splittings between the quasi-degenerate states, while the off-diagonal elements proportional to e 
contribute only to the active-plus-sterile mixing matrix. 
The 6x6 neutrino mixing matrix is 



V 





Ve2 






K2' 




\ 






^,3 










Vr2 


VrZ 


VrV 


Vr2' 


Vry 












Vs,2' 












Vs^V 


Vs,2- 




I 




Vs,2 


^S33 




Vs,2' 





^ / 13 -b 



u -u 
13 13 



13-V 

l3 + <^- 



(11.9) 



(11.10) 



where we define the weak eigenstates as Va, ct — e, fi, r, si, S2, S3 and the mass eigenstates as i/^, i = 1, 2, 3, 1', 2', 3'. In 
the Dirac limit £,„ 0, Vai = —Vai' = Uai/V2 (a = e, ^, r, i ~ 1, 2, 3) and rrii = —m[ [i = 1, 2, 3). We will order our 
states in the "usual" way [l^: ml < < mf, m^, < rn^i < m^, (m| < < m|, my < m\, < m\,) in the case of a 
normal (inverted) mass hierarchy. We assume the parameters small enough that the same mass hierarchy applies 
for the primed and unprimed eigenmasses. More specifically, |mf — m^, | ^ mf for all i = 1,2, 3. 

In the next section we will concentrate on two subsets of the six neutrino mixing scenario: the case of two active 
and one right-handed neutrino, and the case of two active and two right-handed neutrinos. These will be described 
in more detail below and should serve as more concrete pedagogical examples. 

In the case of two active {u,, and j/q, a linear combination of and Vt) and one sterile neutrino Ug, the 3x3 
neutrino mass matrix can be written as 
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(11.11) 



where Cm = me is the Majorana mass of the right-handed neutrino. In the limit e ^ 1 (and ignoring the case when 
is very close to or 7r/2), and assuming all parameters are real. 
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This system is described by three neutrino mass eigenstates: a massless one, vi, which is a linear combination of and 



and two massive ones, V2^2'i which are almost degenerate in mass-squared in the limit e <C 1: 
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m| — m^i = 2m?e. In vacuum, the survival probability of electron neutrinos P^e with energy E after a distance L has 
been traversed is 
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On the other hand, the oscillation probabilities of active neutrinos into sterile neutrinos are 

'2m^eL'" 



Pes = sin^ 6 sin^ 
Pas = cos^ 6 sin^ 
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Eq. (jII.14p agrees, of course, with the well-known two- neutrino oscillation probability in vacuum in the limit e ^ 1. 
On the other hand, in the very long distance limit, 2m?eL ^ E, the oscillation due to the small mass-squared splitting 
can be observed on top of the averaged-out "active-active" oscillations. 
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In the case of two active {ve and a linear combination of and Vr) and two sterile neutrinos Vs^ and Vg^, 
4x4 neutrino mass matrix can be expressed as follows, assuming all parameters real: 
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For concreteness, we allow 6 G [0, 7r/2] and define m| > mj. In this case, re-expressing tmi = wiei, em2 = fnif-i, 
£m3 = [("^2 ~ nT-i)/nT-2]^3, and assuming all ^ 1 and 9 not too close to or 7r/2, 
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where 
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and the 2x2 matrices 
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This system is characterized by four mass eigenstates vi^i'2, i^i'i 1^2' which are pair- wise quasi-degenerate in mass- 



squared. Oscillations are described by six distinct oscillation frequencies: Ami2 = ^2 — 
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2m|e2 (in the limit 61,62 <C 1). It is convenient to also express the two small frequencies in terms of the 
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elements of the Majorana mass matrix for the right-handed neutrinos: Am^^ = 2?7iie„ii, and Ar 
all elements of 6^ are of the same order of magnitude, one expects ATO2/2 > Am^-^ since m^ > mf. 
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In vacuum, the survival probability of electron neutrinos P^e with energy E after a distance L has been traversed is 
sin^ d ej\ . o f /\m'^oL\ f ^ €i 
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(11.23) 

: 1,2,3 (ef,eie2; etc). It is trivial to note 
that, for small enough ei,e2,e3 the two-flavor vacuum oscillation expression is reproduced, as expected. Finally, in 
the very long baseline limit, assuming that the fast 12 (and 1'2, 12', 1'2') oscillations average out, 
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Above, O(e^) indicates terms which are of order the product of two e; 
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Again as expected, in the very long-baseline limit the electron neutrino survival probability is equal to the averaged 
out "active oscillation" effect, plus two long- wavelength components driven by the I'l and 2'2 mass-squared splittings. 
Two features are readily visible: one is that there no dependency on the off-diagonal es parameter (this only appears 
as O(e^) level corrections to the coefficients of the different terms in Eq. (|TL24|). The second is that, assuming 
cos^ 6 ^ sin'* 6, we expect to be more sensitive to £2 than ei, as the 2'2-oscillations "turn on" before the 1'1-oscillations 
(remember 7712 > ra\). 



III. CURRENT CONSTRAINTS 



If one assumes that there are three right-handed neutrinos with Majorana masses much smaller than the Dirac 
neutrino masses that govern "active" neutrino oscillations, current experimental data can be used to constrain the 
right-handed neutrino Majorana mass matrix £„. There are, however, too many free parameters in and the 
individual upper bound on each of these is not very illuminating and outside the aspirations of this work. Such a 
bound will also depend on aspects of neutrino masses and mixing that are currently unknown, including the neutrino 
mass hierarchy and |J7e3p- We are, however, interested in asking what is the constraint on assuming that all its 
elements are of the same order of magnitude. Hence we will concentrate on the current experimental upper bound on 
the element of tm that is best constrained. This will become clear in the following paragraphs. 

If the three neutrino mass eigenstates 1^1,1/2,7^3 identified experimentally [l^ are in reality "split" into, say,* six 
states 7^1, j/2, 1^3, 1^1', 1^2', ^'3'- As discussed in the previous section, these splittings will manifest themselves via new, 
very long wavelength oscillations characterized by the A77i^,j mass-squared differences (i = 1,2,3). These, on the 
other hand, are proportional to emaini and the largest new mass-squared splitting is associated to the e,„.. value 
associated to the largest 777^. In turn, the largest 777^ value depends on the neutrino mass hierarchy. In the case of an 
inverted mass hierarchy, m\ > m\ 2 x 10~^ eV^ is the largest 777^. On the other hand, if the hierarchy is normal, 
777! > 2 X 10^'^ eV^ is the largest 777? while 7773 > 777^ and 7712 ^ 8 x 10^^ eV^. 

Ignoring |t/e3|-driven effects, A7773,3 is best constrained by the disappearance of muon-type neutrinos and antineu- 
trinos produced in the atmosphere. We estimate that these experiments are sensitive to Cm values that lead to new 
oscillation lengths which are not much larger than the diameter of the earth, or 

L^2 - - f lE!^^ ( ^ ^ X 102 km < 10^ km. (III.l) 

£3 V "^3 /V 100 MeV / 



* The current data are also consistent with two right-handed neutrinos. In this case, the massless neutrino mass eigenstate is not split, 
and also does not contain a sterile neutrino component. This is identical to the 2+1 case [Eq. IIII.lll l] discussed in the previous section. 
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Hence atmospheric experiment s can "see" £3 values larger than around 10~^. This, in turn, translates into a sensitivity 
to Ej^^^elements of order ^ vlO^ x 10~^ eV. A more detailed estimate can be extract from the analyses performed 
in [201 ■ In more detail, atmospheric data constrain the new mass-squared difference to be less than about 10~^ eV^. 
In the case of a normal mass hierarchy, this translates into Ae^ x 10~^ eV^ ^ 10~^ eV^ or £3 ^ 0.03 and £„ ^ 0.001 eV. 
The estimate above was made for m| = 2 x 10~'^ eV^. 

Am2/2 (and Am^,j^) is best constrained by solar neutrino experiments. Ultimately, one is sensitive to oscillation 
lengths of order the earth-sun distance (1 A.U.= 149.6 x 10^ km). Naively, the sensitivity to £2 and ei can be 
estimated from 

which translates into £1^2 ^ 10^^ if both mi and TO2 are of the same order. This translates into a sensitivity to 
£to elements of order £„ > V8 x 10^^ x lO"*" eV. Both sensitivity estimates (£3 ^ 10"^ and £1^2 ^ 10~^) are 
quoted assuming the Dirac neutrino mass hierarchy is normal. In the case of an inverted hierarchy the sensitivity 
of atmospheric neutrino experiments to £3 is markedly worse, while that of solar neutrino experiments to £1,2 is 
markedly better. Finally, if the Dirac neutrino masses are quasi-degenerate and much larger in magnitude than 
^yArnf^ ~ 0.05 eV (the atmospheric mass-squared difference), the sensitivity to all ei is expected to be better than 
estimated above. 

In conclusion, constraints on £1,2 from solar neutrino data provide, by far, the best bound on a given e^n element, 
regardless of the neutrino mass hierarchy. Therefore, we will concentrate on the effect of the seesaw right-handed 
neutrinos in solar oscillations. We will further assume {Uesl is small enough so that 1/3 and related effects in 
experiments with electron neutrinos in the initial state are negligible and will concentrate on an effective system 
consisting of the electron neutrino, the linear combination Va of Vp, and v^- orthogonal to the one in ^'3 3/ and either 
one or two sterile states. 



A. Oscillation of Solar Neutrinos, 2-\-l case 

We first discuss the case of only one right-handed neutrino, Eq. (jll.lip . where toi = 0. The 2-1-1 scenario captures 
most of the physics of the 2-|-2 case, which will be discussed in the next subsection, and contains all the relevant 
information in the limit mi <C m2. In the 2 -|- 1 case, the effect of right-handed neutrinos is entirely captured by 
one dimensionless parameter e. Furthermore, the neutrino mass responsible for the oscillation of solar neutrinos and 
reactor antineutrinos at KamLAND is uniquely determined: m^ = (7.59 ± 0.21) x 10~^ eV^ [2l| from KamLAND 
data. We will show that e is constrained to be small enough that this measurement is not affected by the presence of 
the sterile neutrino (sec Eq. pi.l4[) '). The other parameter in M^, the mixing angle 6, is mostly determined by solar 
data, more specifically those from SNO and Super-Kamiokande. We will argue later that the impact of £„ on the 
determination of 9 is not significant. 

We first briefly discuss how neutrinos produced in the core of the sun propagate towards the surface of the sun 
and then to a detector on earth, and calculate the relevant transition probabilities. We hence construct the effective 
Hamiltonian H in the flavor basis in the presence of the matter [2^ : 




m sm y cos u me sm t 



2E 2E ^2E 

sin 6 cos ^ m cos m"ecos6 
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2E 2E 2 

rn^esinO rn^n cos m (1+e ) 

2E 2E 2E 



(111.3) 



where A = v^Gi^A^e is the matter potential due to electrons, Gp is the Fermi constant, is the position-dependent 
electron number density, and E is the neutrino energy. In principle, we should also consider the effect of neutral 
current interactions due to the presence of neutrons. This contribution, however, is expected to be negligible in the 
sun given that it consists mostly of hydrogen (hence "neutron-poor"), and the proton and electron neutral current 
contributions cancel out. Under this approximation, analytical solutions for the transition probability can be easily 
obtained as long as Ng is a simple enough function of the neutrino position. 

It is easy to understand how an electron neutrino produced in the sun's core propagates first to the sun's surface 
and later to the detectors on earth. It is illustrative to discuss what happens in the £ ^ limit. In this case, one 
of the eigenvectors of H (Eq. pil.3[) ). with eigenvalue m? /2E, is independent of A and is purely sterile. In the 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

R(in units of the radius of tlie sun) 



FIG. 1: The evolution of instantaneous eigenvalues of the Hamiltonian inside the sun in the 2+1 scenario, Eq. (|III.3[I . for 
E = 5 MeV, = 8.1 x 10"'' eV^ sin^ 9 = 0.3, e = 3 x lO'V 



basis defined by Eq. (jII.13p it is neither 1^2 nor 1^2' but l/v^(i^2 + 1^2')-^ The other two states are purely active and 
can be obtained by diagonalizing the familiar matter-affected two-by-two neutrino propagation Hamiltonian. For 
values of interest, electron neutrinos arc produced in the sun's core and propagate adiabatically, exiting the sun as an 
incoherent mixture of vi and l/-\/2(i^2 — 1^2' ) with probabilities cos^ Bm and sin^6'M, respectively, where the matter 
mixing angle is defined at the production region by the familiar expression 

• Asin26' „^ Acos26'-A 

sm2gM = ^ cos2(^M = ^ (III.4) 

' (A cos 261 - A)2 + A2 sin^ 29 J (A cos 26 - A)^ + A^ sin^ 261 



where A = /2E. In this case, P^e at the surface of the sun (or the earth) is given by 

lim Plf^' ^ cos" 0m|14i I' + \V,2 - Ve2' I' = cos" Om cos" e + sin" Om sin" 9, (III.5) 

e— >0 2 

which is the familiar result for solar neutrino oscillations in the so-called LMA region. 

Figure [T] depicts the evolution of the instantaneous eigenvalues of H inside the sun, as a function _R, the distance 
from the sun's center, for e 7^ 0. For illustrative purposes wc choose, for the Hamiltonian parameters, E = 5 MeV, 
to" = 8.1 X 10~^ eV", sin" 9 = 0.3, e = 3 x 10^^ The LMA MSW resonance (between the largest eigenvalue and the 
smallest one) can be readily identified (at R ~ 0.1). One can also see the impact of e 7^ 0: the top two eigenvalues are 
"split" for all values of R. Furthermore, there is another "resonance" [2^ due to the small but non-zero e at larger R 
between the two heaviest Hamiltonian eigenstates. This will be discussed momentarily. 

Assuming that the crossing of the LMA MSW resonance is adiabatic and that e ^ 1 , the electron neutrino survival 



t In the limit e — > 0, 1^2 and U2' have the same mass-squared and hence any Hnear combination of them is an cigenstate of the propagation 
Hamiltonian in vacuum. 
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probability can be written as 
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(III.6) 



1 - e-T 



7 ■ 



9.8 



(lO-O 



8 X 10-5 eV2 



0.862 MeV\ 
E 



(111.7) 



V is given by Eq. pi.l3p and Vmat is the unitary matrix that diagonahzes Eq. pil.3p at the production point. Eq. pil.6[) 
can be understood as follows. The electron neutrino is first expressed in the basis of the Hamiltonian at the production 
point. It then propagates adiabatically (each component acquiring a phase factor c/ii, i = 1,2,3) until it gets to the 
location of the resonance between the two "heaviest" instantaneous Hamiltonian eigenstates. Around that point, the 
unitary evolution is characterized by the matrix containing which is the crossing probability between the two 
Hamiltonian eigenstates. Finally, it propagates from that point until it gets to the earth (each component acquiring 
a phase factor cfi'^, i = 1, 2, 3), where we compute the probability that this state is an electron neutrino. For a similar 
discussion of Pee in the presence of a new small mass-squared difference, see, for example, [20| . 

We can estimate Pee for e ^ 10"'^ when Pc in Eq. pil.7p vanishes to a very good approximation for all neutrino 
energies below 10 MeV. If we consider the case ^ ^ A in the production region (this is an excellent approximation in 
the upper energy range of the solar neutrino spectrum). Pee is very easy to compute since the electron neutrino is 
a Hamiltonian eigenstate at birth (corresponding to the "heaviest" Hamiltonian cigcnstate) and the entire evolution 
of the state inside the sun is adiabatic. In summary, the electron neutrino exits the sun as the heaviest neutrino mass 
eigenstate, which we will assume to be V2 (which is the case for e positive) and 
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It is easy to see that under these circumstances there is no value of 9 that provides a good fit to the solar neutrino 
data, as we will discuss in more detail in Sec. IIII CI 

On the other hand, in the limit e ^ lO""^, it is sufficient to keep only the leading order (e 0) limit of Vmat- In 
more detail, 



VI 




(III.ll) 



where 6m is defined in Eq. pil.4p . At the surface of the sun, after integrating over the neutrino production region, the 
electron neutrino is an incoherent mixture of Ui (with probability cos^ 9m) and y/T' 



Pc'^2 — \fP'cU2''^ (with probability 
sm^QM)- The latter may undergo vacuum oscillations on its way between the sun's surface and the earth with an 
oscillation probability given by 



P(V/1 - PcV^ - ^JPcVT ^ Ve) = 
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sm 



2m^eL 



l + 2v/Pe(l-Pe)cOS 



2E 

2m'^eL 
2E 



(IH.12) 



t An irrelevant relative phase factor between V2 and V2' been omitted. 
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FIG. 
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2: PeC for solar neutrinos at the surface of the earth, — e,a, s, as a function of neutrino energy for — 7.6 x 10 ^ eV^ 
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Combining all the information above, 



Pee ~ cos^ 6m cos^ 9 + sin^ 6m 



Pea ~ cos^ 6i\i sin^ I 



Pe, 



sm Um-^ 



sm Um - 



sm f 
2 

cos^ 6 



1-2/^(1^^ cos 



l + 2v/Pc(l-Pc 



/2m2eL 

V 2^; 



l + 2v/Pc(l-Pc)cos 



2m2eL 

2^; 

2E 



(III.13) 



Figurc[2]dcpicts Pg^ C = e, a, s for solar neutrinos as a function of energy for m? = 7.6 x 10~^ eV^, sin^ 6 = 0.31 and two 
different values of e = 1 x 10~^ and e = 5 x 10~®. In order to compute PeC; we integrate over the neutrino production 
region in the sun's core, taking into account that 6m depends on Ne at the production point. For high energy solar 
neutrinos Pes is small and Pee is very similar to the standard LMA solution to the solar neutrino puzzle. For low 
energy solar neutrinos, the oscillatory pattern becomes more pronounced and Pes becomes significantly nonzero. As 
already briefly discussed, for much larger values of e, Pc ^ and Pes — sit? 6m 1^2 is roughly equal to one half for 
high energy solar neutrinos. On the flip side, for smaller values of e, Pc 1/2 and the argument of the cosine in 
Eqs. (jlll.lSP is vanishingly small. In this case, as already discussed, sterile neutrino effects disappear and we recover 
the standard LMA solution to the solar neutrino puzzle. 



B. Oscillation of Solar Neutrinos, 2+2 case 

It is straight forward to generalize the results obtained in the previous subsection to the case two active (ve and 
Va) and two sterile neutrinos (y\ and V2). In this case, the Hamiltonian that governs neutrino oscillations in neutral 
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FIG. 3: The evolution of instantaneous eigenvalues of the Hamiltonian inside the sun in the 2+2 scenario, Eq. (|III.14|) . for 
= 10 MeV, m? = 1 X 10"® eV^ mi = 8.2 x 10"^ eV^ sin^ = 0.32, ei = 3 x lO'S £2 = 4 x 10"\ and e„3 = 10"^ eV. 



hydrogen matter characterized by an electron number density N,, is 



H = 



2E 



A 



Acos^ 
A cos sin 9 



2E 



A cos 9 sin 9 
Asin^^ ' 



2E 



V 



2E 



2E 



63 Am~ 



2E 
2B 



0(6^) 

m^(l+0(e")) 
2_B 



Here, My is given by Eq. (jll.lSp while U2 and I2 are defined in Eq. pi.2ip . We define Ar 



(III.14) 

i\. Note that a 

diag(J72, 12) rotation expresses the Hamiltonian in an active-sterile mass-squared basis where the top two components 
consist of two mass-squared states that are mostly active, while the bottom two components consist of two mass- 
squared states that are mostly sterile. The Hamiltonian is block diagonal in this basis in the limit £1,2,3 0- In this 
limit, of course, we are left with the standard MSW Hamiltonian for two active neutrino flavors. 

Figure [3] depicts the evolution of the instantaneous eigenvalues of H inside the sun, as a function R, the distance 



from the sun's center, 



£1,2,3 7^ 0. 
.2x10-5 eV^, sin^e = 0.32, £1 = 



for 

.2 , 



For illustrative purposes, we choose £^ = 10 MeV, 



3 x 10-1, £2 = 4x 10-\ and e 



= 1 x 10-^ eV'', 
= 10-'^ eV for the Hamiltonian parameters. 
The LMA MSW resonance (between the largest eigenvalue and the smallest one) can be readily identified (at i? ~ 0.1). 
One can also see the impact of £1,2,3 7^ 0: the top two and the bottom two eigenvalues are "split" for all values of R. 
Similar to the 2-1-1 case discussed in the previous subsection, there are two other resonances [23 due to the small but 
non-zero £1,2,3 at larger R between the two largest and two smallest Hamiltonian eigenstates. From the Hamiltonian 
Eq. (jIII.14p . it is easy to see that these are governed by the small parameters £1 (resonance between the two lightest 
states) and £2 (resonance between the two heaviest states). £3 connects the heaviest "active" state to the lightest 
"sterile" state (and vice- versa) and, for physically interesting values of the neutrino oscillation parameters, its effects 
are not visible. 

As in the 2-1-1 case, Am^ = (7.59 ±0.21) x IQ-"^ eV^ is fixed by the results of the KamLAND experiment, and sin^ 9 
is constrained to be large. In this case, electron neutrino production in the sun's core and subsequent propagation can 
be described as follows. The electron is produced as an incoherent superposition of the two mostly active states (this 
is true as long as A ^ £1,3^1/ {2E) for all values of E of interest). The two components then evolve adiabatically until 
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they hit the resonances governed by ei and 62. Hence the electron neutrino exits the sun as an incoherent mixture of 
Vl — Pcii^i — ^Pc\vv (with probabihty cos^ 6*^/) and \J\ — PciV2 — \JPc2V2' (with probabihty sin^ 6*71/)-^ Here, Qm is 
defined by Eq. pff.4p . The crossing probabihties are given by 



Pr. 



1 -e-T 



10- 



X 10-5 eV2 



0.862 MeV\ 
E ) 



1,2, 



(HI.15) 



where Vg^i are given in Eq. pi.20p . On their way to the surface of the earth, both of these states can undergo vacuum 
oscillations driven by the mass-squared differences 2m\t\ and 2m^e2 [see Eq. (|II.24|) ]. One can write 



Pee -C0S2 X P{^J\ - Pd 



'P'lyv V,) + sin^ X P(v/1 - Pc 



2V2 



'P^2V2' -*Ve), 



(IH.16) 



where P{\/\ — Pc\V\ — \/Pc\Vv 



i^e) is given by Eq. piI.12p with Pc, m, e replaced by Pci, nT-i, ei, respectively, and 
Ue2,e2' replaced by Uei,ei' , respectively. Similarly, P(-\/l — Pc2i^2 — \JPc2V2' v^) is given by Eq. piI.12p with Pc,m,e 
replaced by Pc2, 1^2,^2, respectively. After the dust settles, 



cos 



l + 2v/Pci(l-Pci)cos 



COS 6m- 



sm 



l + 2v/Pci(l -Pci)cos 



E 
•m\eiL 



sin^ 9m - 



sm 



l+2v/Pc2(l-Pc2)cOS 



+ sin^ 9m- 



cos 



r— — -, — — 7 / mie^L 

1 - 2v/Pci(l-Pci)cos f 



SIT? 9m-z 



l + 2v/Pc2(l-Pc2)cOS 



m%e2L 



1712^2 L 

E 
m%e2L 



E 



1 - 2v/Pc2(l-Pc2)cos 



E 



(IH.17) 



Pee = COS 
P 

^ ea 

Pes = COS^ 9m ^ 

where Pes = Pe,5i + Pesa and 9m is given by Eq. pll.4p . In all expressions above, we have assumed all 61,2,3 much 
smaller than 1. This will be justified a posteriori (see Sec. IIII C|) for ei and £2 and we postulate that it is valid for £3. 
The impact of a "large" £3 is mostly felt in the diagonalization of H [Eq. piI.14p ] in the production region, which we 
assumed is properly described only by the matter mixing angle 9m [see Eq. pil.llj) ]. The bound one would obtain on 
£m3 via these effects is much worse (by many orders of magnitude) than the bounds on £,„i and 6^2 discussed in the 
next subsection. As explained earlier, we are interested in the most stringent bound on any element of £m and will 
not consider these £3 effects henceforth. On a related note, Eqs. piI.17|) do not depend on £3, similar to Eq. pi.24[) . 
For this reason we will not discuss £3 effects or bounds. 

Figure U] depicts Pe^ for solar neutrinos at the surface of the earth, C = 

Qj^ S ^ clS Oil function of neutrino energy for 



1.0 X IQ-^ eV^ 



5.2x10"'^ eV^ sin^ 



0.31 and £1 = 1 x 10 °, plus two different values of £2 = 1 x 10 



and £2 = 1 X lO"''. It is easy to see that dominant £1 effects are most visible for low energy solar neutrinos, when 
cos^ 9m is significantly different from zero [Fig. [4l^T0P)]. On the other hand, dominant £2 effects are identical to 
the ones observed in the 2+1 case, where £2 plays the role of e. The main distinction is the fact that in the 2+1 
case was constrained by KamLAND, while in the 2+2 case neither ml nor are strongly constrained by data, 
except through their difference: ra^ — m\ = 7.6 x 10"^ eV^. In the case when both £1 and £2 are non-negligible 
[Fig-SlBOTTOM)], one can see that £2 effects are dominant for high energy solar neutrinos (c/. Fig. [2]), while £1 and 
£2 effects "add up" for the low energy solar neutrinos. 



C. Constraints from Experimental Data 

The flux of electron neutrinos from the sun has been measured, sometimes as a function of energy, by a variety 
of experiments and in a variety of ways. The flux of active neutrinos from the sun has also been unambiguously 
measured for *B neutrinos. Here, we use recent data from [2^ [2^, [l^, [H, [2^, [sO, [lH in order to constrain e (in the 
2+1 case) and £1,2 (in the 2+2 case). In order to compute the expected event rates at the different experiments we 
use parameters from the BS05(OP) version of the Standard Solar Model (SSM) [3^ . 

We do not attempt an exhaustive analysis of all neutrino data but instead concentrate on the features that most 
impact our results. We make use of the following observables. More details are provided in the Appendix. 

• (j>^'^{SNO), the integrate ®B neutrino flux measured via J^e+^H e^+p+p at SNO [H, [H, [131 . This observable 
is sensitive only to the electron neutrino component of the solar neutrino flux. 



5 Relative phases between Ui and i/j/, j = 1, 2 have been ignored as neither of them is observable in practice. 
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FIG. 4: PeC for solar neutrinos at the surface of the earth, ( = e, function of neutrino energy for ml — 1.0 X 10"'' eV^ 



mi = 8.2 X 10" 



eV^ sin^ 6 = 0.31. TOP: ei = 1 x 10"^ £2 = 1 x 10"^ BOTTOM: ei = 1 x 10"^ e2 = 1 x 10~ 



(^^'-^(SNO), the integrated neutrino flux measured via i^e.a+^H ^ ^^e.a + n + p at SNO 
observable is sensitive to the active neutrino component of the solar neutrino flux. 



mm 



This 



lES 



(SuperK), the *B neutrino flux measured via 



at SuperKamiokande [28[ integrated 



over a handful of recoil electron energy bins. This observable is sensitive to the active neutrino component of 
the solar neutrino flux. The cross section for muon or tau-type neutrino-electron scattering is about 0.16 times 
that of electron-type neutrino-electron scattering when integrated over the neutrino energy spectrum. SNO 
also has data on neutrino-electron elastic scattering (also included), but with significantly less statistics than 
SuperKamiokande. 

• (^^"^(Borexino), the ^Be neutrino flux^ measured via i^e.a + e" i^e,a + at Borexino [3l|. This observable 
is sensitive to the active neutrino component of the solar neutrino flux. The cross section for muon or tau- 
type neutrino-electron scattering is about 0.21 times that of electron-type neutrino-electron scattering for ^Be 
neutrinos. 

• 4>''^°', the measurement of the solar neutrino flux using inverse /3-decay in Gallium [29l |. This observable is 
sensitive to neutrinos produced in all distinct neutrino-producing fusion reactions, but only to the electron 
neutrino component of the solar neutrino flux. 

• (/)*-^', the measurement of the solar neutrino flux using inverse /3-decay in Chlorine [30| . This observable is 
sensitive to neutrinos produced in all distinct neutrino-producing fusion reactions, except those involving the 
p + p reaction, but only to the electron neutrino component of the solar neutrino flux. 

All flux measurements are defined assuming that the incoming neutrinos are all electron neutrinos. For example. 



lES 



(Borexino) 



P.. 



- P 



(III.18) 



^ Borexino is also sensitive to neutrinos produced in the CNO cycle. These account for around 10% of the neutrinos with energies around 
those of the ^Be neutrinos and have been neglected. The current Borexino data is consistent with a vanishing CNO neutrino flux |31|| . 
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where <j>7Bc is the SSM expectation for the ^Be solar neutrino flux and cr^^e is the cross section for elastic C-type 
neutrino-electron scattering, ^ = e, a. 

Before presenting the results of a fit to the data spelled out above, it is illustrative to describe the dominant 
aspects of the data that constrain the presence of the extra mass-squared differences. For ^B neutrinos, the combined 
SNO and Super-Kamiokande data not only reveal that P^e ^ 0.3 for neutrino energies above a few MeV, but also 
reveal that Pea ^ 0.7. A more detailed analysis, presented in the Appendix, reveals that for ^B neutrinos. Pes < 0.37 
at the 3(7 confidence level. As advertised, this rules out at least e > 10"'^ in the 2+1 case [see Eq. pil.lOp ]. Upon 
further scrutiny, it is easy to see that, for e > 10~^, Pes ~ 0.5 in the 2+1 case, which is ruled out by the current data. 

For solar neutrino energies below 1 MeV or so, information regarding neutrino oscillations is dominated by Borexino, 
the Homestake experiment, and the Gallium experiments. The Borexino experiment measures, with good precision, 
mostly for '^Be solar neutrinos {E^, = 0.862 MeV), Pee = 0.56 ± 0.10 in the limit Pes [IH. For larger values of 
Pes, Pee is driven towards higher values in order to compensate for the depletion of the Ua component of the ^Be solar 
neutrino flux at the earth. This effect is depicted in more detail in the Appendix. For lower energy neutrinos, the 
Gallium experiments, taking into account the data, mostly, from Borexino, SNO, and SuperKamiokande, constrain 
Pee ~ 0.55 for pp neutrinos, and are insensitive to whether those convert into active or sterile neutrinos.** 

Figure O depicts Ax^ = ~ X^i^ = 0) as a function of e, in the case of two active and one sterile neutrino, for 
= 7.6 X 10^^ eV^ and sin^ 6 ~ 0.31. is the result of a x^ fit to the data described above, and we will use it in 
order to establish an upper bound on e. More specifically, we will state that values of e associated to Ax^ > 4, 9 are 
ruled out at the two, three sigma level. While we quote a bound for a fixed value of , sin^ 9, we have verified that 
a very similar bound is obtained for different values of these oscillation parameters (in which case we also include 
the measurement of from the KamLAND experiment in the x^ function). We note, for example, that there is no 
good fit for values of 9 outside the currently best fit region for sin^ 9i2 [2l| , even if one considers significant values of 
e. From the figure, we estimate that the largest value of e allowed at the three sigma level is e < 2.0 x 10~^, while 
at the two sigma level e < 1.2 x 10"''. Furthermore, the region 2.4 x 10~^ < e < 7.1 x 10~^ is ruled out at the two 
sigma level. In this region Pes is large for solar neutrino energies around 1 MeV, while Pee is less than one half, in 
contradiction with the Borexino data. For slightly larger values of e, the peak in Pes moves to solar neutrino energies 
larger than 1 MeV but smaller than the threshold for the SNO and SuperKamiokande experiments, where there is 
virtually no experimental information. In this case, the fit to the data is just as good as the one obtained for e = 0. 

It is also worthwhile pointing out that, for e ~ 1 x 10~^, Pee for ^B neutrinos is expected to be mostly constant 
as a function of energy, contrary to the standard LMA case, where a slight decrease of Pee as a function of energy is 
expected. This behavior is observed for the appropriate parameter choices in other "new physics" scenarios including 
non-standard neutrino interactions [11] , mass- varying neutrinos [l^ , and the introduction of light sterile neutrinos 
that mix very weakly with active ones (35j . Such a behavior was advertised as solution to the fact that the ^Be 
electron neutrino flux, extracted using the Chlorine data, was lower than the one predicted by the canonical LMA 
solution to the solar neutrino puzzle. This tension in the data has been significantly relaxed with the introduction of 
those from Borexino. 

The bounds above translate into 2TO^e < (1.8, 3.0) x 10~^^ eV^ as the two, three sigma upper bound on the induced 
small mass-squared difference and e™ < (1.0,1.7) x 10~^ eV at the two, three sigma level. Note that our bounds 
agree qualitatively with those obtained, under slightly different circumstances and after some re-interpretation of the 
relevant observables, in (2^. The analysis in did not include results from Borexino, not available before 2008. 

In the 2+2 case, we wish to place bounds on both ei and £2, and discuss which between emi = wi^i and 6^2 = 'm2e2 
is more severely constrained. Unlike the 2+1 case, the data (mostly from KamLAND) do not specify mf or m|, but 
constrain m| — = 7.6 x 10~^ eV^. As discussed in the previous subsection, the solar data constrain mfei and 
m\e2- Keeping all this in mind, we will first discuss the bounds of ei and £2 for fixed m\ = 1.0 x 10~^ eV^ and 
m\ — 7.7 X 10~^ eV^ and then discuss the bounds on for various values of m\^m\ by varying these two masses 
while maintaining m| — constant and in agreement with the KamLAND data. As in the 2+1 case, we will discuss 
bounds for fixed m\ — rnf — 7.6 x 10"^ eV^ and sin^ 9 = 0.31. We have checked that a very similar bound is obtained 
when one allows these two parameters to vary in the fit. 

In more detail, we discuss the bound on ei in the limit £2 very small and vice- versa. The reason for doing this is that 
£1 X £2 correlated effects are not very large. In the limit where £1 effects are negligible, £2 bounds are similar to those 
on £ in the 2+1 case. In more detail, 2m2£2 < 3.0 x 10~^^ eV^ at the three sigma level. In the limit where £2 effects 
are negligible, we obtain, following a procedure identical to the one discussed in the 2+1 case above, £1 < 5.3 x 10~^ 
at the two sigma level and £1 < 8.8 x 10^^ at the three sigma level for mj = 1.0 x 10~^ eV^. This translates into a 



* The bound on P^e depends indirectly on Pes for higher energy neutrinos as these modify the measurement of Pee for the high energy 
neutrinos and hence the extraction of Pee for pp neutrinos from the Gallium data. Details arc provided in the Appendix. 
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FIG. 5; - X^(e = 0) as a function of e for = 7.6 x 10"^ eV^ sin^ 6 = 0.31, in the 2+1 case. 



new mass-squared difference 2m\ei < (1.1, 1.8) x 10 eV^ at the two, three sigma level. 



We are interested in upper bounds for the elements of 
we obtain, at the three sigma level, 



From the bounds on the mass-squared differences above. 



e„i < 8.8 X 10-10 cV (i^) 
< 1.7 X 10-9 evfM><ioii^ 



(III.19) 
(III.20) 



These bounds are obtained for sin 9 = 0.31 and — ml = 7.6 x 10^^ eV^ and are hence correlated. In the case of an 
inverted mass hierarchy, for example, if mi = 50.000 x IQ-"^ cV while m2 = 50.076 x 10-'^ eV, the three sigma bounds 
above translate into Cmi < 1-8 x 10"" eV and 6^2 < 3.0 x 10-^° eV. On the flip side, while m^ > 7.6 x 10-^ eV^ 
is bounded from below, mi can be arbitrarily small, which indicates that the upper bound on e,„i, while naively 
stronger than that on 6^2, could be significantly worse. 

Before moving on to current non-solar and future constraints on Cm, wc would like to emphasize a few relevant points. 
We have not considered earth matter effects while computing our oscillation probabilities, neither have wc worried 
about large seasonal variations of the solar neutrino flux induced by the very large oscillation length induced by the 
new, very small mass-squared differences. The fact that we can ignore day-night effects is justifled a posteriori. Earth 
matter effects are most significant when the matter potential of the earth is of order Am'^ /2E. For ^Be neutrinos 
{E = 0.862 MeV) "large" (more than one percent) day-night effects are expected for Am^ > 10^^ eV^ (see, for 
example, [11]) so we estimate that, for solar neutrino energies above 100 keV (below the threshold of the Gallium 



Strictly speaking, some of the approximations made in Sec. |ll] fail when mi is of order emi and hence the bound Eq. piI.19| l is not 
applicable when mi < 10~® eV. In this case however, the generic upper bound on em, as discussed in the introduction to Sec. IIIII is 
dominated by Eq. Jill. 201 1. For even smaller values of mi < 10~^ eV we anticipate that there are no relevant effects due to the fact that 
our approximations fail. On the contrary, in this case we expect oscillations due to sterile effects to be entirely dominated by the 2+1 
case discussed here. 
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experiments) large day-night effects are expected for Am^ ^ 10~^ eV^. Tfiis is two orders of magnitude larger 
than the worse upper obtained for jml ~ "^2' I- These estimates are in agreement with more detailed computations 
performed in Seasonal variations, the other hand, are expected for < 1 MeV if the bounds on are saturated. 
These could be observed in the at, for example, the Borexino experiment. For an estimate of how large an effect one 
may hope to observe, we point readers to [2(]|.[37|. 



Whether the neutrinos are Dirac or pseudo-Dirac neutrinos is a notoriously difficult issue to resolve, especially in the 
limit when the right-handed neutrino mass matrix e„i is tiny. Other than observing the long-wavelength oscillation 
driven by the new mass-squared differences between the quasi-degenerate mass-squared states, it seems virtually 
impossible to construct observables that will reveal the pseudo-Dirac nature of the neutrinos. 

It is well-known that for pseudo-Dirac neutrinos which arise from a seesaw Lagrangian [Eq. (jl.l[) ] the neutrino 
exchange contribution to neutrinoless double-beta decay vanishes almost perfectly [T^ . One way of understanding 
this is to note that the contributions of the different quasi-degenerate states to neutrinoless double-beta decay are 
equal and opposite, and cancel pairwise. Another way of appreciating this fact is to see that, when all neutrinos are 
lighter than tens of MeV, the rate for neutrinoless double-beta decay is proportional to the ee-element of the full 
neutrino mass matrix, which, as one can quickly read off Eq. (jll.ip . is exactly zero. This "property" is present even 
when the neutrinos are not pseudo-Dirac states, and depends only on all right-handed neutrino masses being small 
enough [1, Q. For the same reason, "all" (for a long list of observables see, for example, [HI) potentially observable 
lepton number violating phenomena are also guaranteed to vanish almost perfectly if the only sources of lepton number 
violation are the Major ana neutrino masses. 

As far as kinematical probes of the neutrino masses are concerned - the most stringent laboratory bounds come 
from tritium beta-decay experiments - the neutrinos behave as if they were Dirac fermions, their "pseudo" nature 
completely obscured by the very stringent constraints from oscillation experiments. The same is true for cosmological 
bounds on the number of neutrino species. For small enough right-handed neutrino masses, the sterile neutrinos are 
too weakly coupled and their states are not significantly populated in the early universe [lo, [H, \^ . 

In the near future, current solar neutrino experiments are expected to improve our understanding of very small 
right-handed neutrino masses. Borexino will accumulate more data (and their systematic uncertainties are expected 
to go down) and improve on their results published in (sij . As mentioned above, Borexino is expected to measure 
the solar neutrino flux as a function of time, and may be sensitive to anomalous seasonal variations. The absence of 
such effects should help erase the local minimum around e = 1 x 10~^ in the 2+1 Ax^ function, depicted in Fig. [5l 
and improve on the bound on e (62) in the 2-1-1 (2-1-2) case. It is also expected that SNO and SuperKamiokande will 
be able to access the low energy end {E < 5 MeV) of the ®B solar neutrino spectrum and shed more light on the 
expected (from the canonical LMA solution) rise in Pee as the neutrino energy decreases. 

Future neutrino experiments are aimed at measuring with more precision (and in real time) the pp solar neutrino 
flux and that of the pep and CNO neutrinos (4l| . These would not only provide a more precise measurement of Pee 
for very low-energy solar neutrinos (the pp neutrinos, E ^ 0.5 MeV) but would also fill the energy "gap" between ""Be 
neutrinos and ^B neutrinos (the pep and CNO neutrinos, i? ~ 1 — 2 MeV). Indeed, Borexino is currently working on 
controlling the cosmogenic ^^C background in order to access the pep and CNO neutrino fluxes. It is expected that 
such results, if consistent with the canonical LMA solution to the solar neutrino puzzle, will improve, perhaps by an 
order of magnitude, the bounds on e in the 2-f 1 case or £1^2 in the 2-1-2 case discussed in the previous section. 

Access to new neutrino oscillation lengths much longer than the earth-sun distance will only be provided by 
astrophysical neutrinos, which travel galactic and extra-galactic distances. In particular, it has been pointed out 
[4^ that studies of the flavor composition of ultra-high energy neutrinos, to which neutrino telescopes are sensitive, 
are sensitive to the new mass-squared differences that characterize pseudo-Dirac neutrinos if these are as small as 
10~^^ eV^ or perhaps 10~^® eV^. Taking these estimates at face value, one would be sensitive to right-handed neutrino 
masses 



where Sm^ is the mass-squared difference which can be potentially accessed by neutrino telescopes p3 | . The sensitivity 
above is estimated in the limit where rn^ = 2 x 10^^ eV^ (in which case the neutrino mass hierarchy is normal) and 
is expected to be more inclusive if all neutrino masses are larger, as was observed in the 2+2 analysis in the previous 
section. 



IV. OTHER CONSEQUENCES AND FUTURE SENSITIVITY 




(IV.l) 
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FIG. 6: Estimate of the magnitude of the mixing between active and sterile neutrinos sin^ 9as as a function of the right-handed 
neutrino mass Mjv, for different values of the mostly active neutrino masses, = 10"\ 10~^ and 10"'^ eV. The hatched 
region qualitatively indicates the values of Mjv that are currently excluded by the world's particle physics data. 



V. SUMMARY, CONCLUDING THOUGHTS 



The addition of gauge singlet Weyl fermions to the Standard Model Lagrangian is a very simple and effective 
way of rendering the neutrinos massive, as required by experimental observations. The most general renormalizable 
Lagrangian consistent with this enlarged particle content, Eq. (jLl[) . while deceptively simple, provides a colorful 
spectrum of different phenomenological consequences, depending on the values of the new parameters in the vSM . 

For any value of the parameters M^- (gauge singlet Majorana mass parameters), it is possible to obtain mostly active 
neutrino masses in the range highlighted by experiments by properly adjusting the value of the Yukawa couplings 
Xai- It is, hence, sufficient to discuss the different phenomenological opportunities for testing the vSM as a function 
of Mij. Further simplifying our discussion, we will consider that all My arc of the same order of magnitude, and will 
refer to all of them as M. 

The largest allowed value of M is of order 10^^ GeV [4^. This comes from the requirement that Eq. (jl.ip is 
appropriate to describe physics at energy scales of order M (unitarity and perturbativity). At these high energies 
(and all the way down to M values of order the weak scale) the only directly observable consequences of Eq. (|I.1|) 
are the small Majorana masses of the mostly active neutrinos. More indirectly, lepton-number violating processes 
involving gauge singlet fermions in the early universe may be responsible for the matter-antimatter asymmetry of the 
universe, through the leptogenesis mechanism [43 |. 

For weak scale M values (1 McV<i M < 100 GeV) the theory contains new neutral heavy leptons, which are, 
predominantly, gauge singlet fermions which should be produced in a variety of weak processes including charged- 
lepton and meson decays (see, for example, [ll|, [l^, [3])- These neutral heavy leptons should also be produced 
in collider experiments and, in principle, their decays should point to lepton number violation at colliders [Zsj . 
Indirectly, the presence of neutral heavy leptons is expected to mediate charged-lepton flavor violating phenomena 
and the apparent violation of universality of the weak coupling constant and of unitarity of the neutrino mixing 
matrix. Whether any of these phenomena is observable in practice depends on the amount of mixing between the 
active neutrinos Va=e,fi,T and the gauge singlet fermions, Vs=si.s2,...- Fig. [S] depicts a naive estimate, as a function of 
M, of this mixing angle sin^ 6as given by Eq. p.ip . In the limit M ^ m^, 

Bin^^as^^, (V.l) 

where m^, arc the mostly active neutrino masses (mi,m2,m3, constrained by oscillation experiments). The figure 
depicts sin^ 9as for three different values of the mostly active neutrino masses, according to guidance provided by 
the current neutrino data. Experiments are sensitive to neutral heavy leptons if siia^ 9as ^ 10~^ (best case scenario), 
which leads one to naively conclude that Eq. (jl.l[) . if responsible for the observed neutrino masses, cannot be tested 
experimentally for M values in this range (1 MeV< M < 100 GeV). It is, however, important to emphasize that 
Eq. (jV.ip can be violated and much larger sin^ 0as values are consistent with the oscillation data [13, HE] ■ Whether 
this possibility is realized in nature can only be established experimentally. 
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Sub-MeV neutral heavy leptons ((1 eV^ M < 1 MeV) can make their presence known in astrophysics, cosmology 
and neutrino oscillation experiments. If M values are in this range, it also turns out that, while neutrinos are Majorana 
fermions, the expected rates for all potentially observable lepton-number violating phenomena is strongly suppressed, 
as the contributions from the mostly active states cancel those from the mostly sterile ones 0, Q . In this mass range, 
mostly sterile neutrinos become an interesting dark matter candidate [l3| and, at the same time, part of this mass 
region is only allowed if one can bypass cosmological constraints on the number of neutrino species 0]. 

In this paper, we were interested in new light neutrinos, or small values of M (zero < Af ^ 1 cV). In this case, 
Eq. (jV.ip does not apply at all, and active-sterile neutrino mixing is naively expected to be large. A/ values of order 
the "active" neutrino masses mi , m2 , are most likely excluded by failed searches for sterile neutrinos in oscillation 
experiments. A detailed discussion of this parameter region is very dependent on the structure of the Yukawa matrix 
A and is not the subject of this paper. Instead, we concentrated in the limit M <C mi, TO2, 

If the gauge-singlet Majorana neutrino masses are much smaller than the Dirac neutrino masses rriai = XaiV, 
where v is the vacuum expectation value of the neutral component of the Higgs doublet, neutrinos are classified as 
pseudo-Dirac neutrinos. When Mij ~ for all neutrinos are Dirac fermions and all experimental data can be 
accommodated perfectly. For Mij ^ but small enough, this conclusion remains true. Our goal was to estimate above 
what small but nonzero value of Mij is Eq. p.ip ruled out by experimental data. We find that the most stringent 
constraints come from the current solar neutrino data. The most conservative upper bound on the best constrained 
element of the gauge singlet Majorana mass matrix, in the basis defined in SecHIl is M < 1.7 x 10~^ eV at the three 
sigma confidence level. This bound is significantly stronger if the neutrino mass hierarchy is inverted or if the lightest 
neutrino mass (toi in the case of a normal mass hierarchy) is larger than 10""^ eV. 

The hatched region in Fig. [S] indicates the values of M for which Eq. (|I.ip is excluded as an explanation for the 
non-zero neutrino masses. It is remarkable that while providing a good quantitative fit to all particle physics data, 
the parameters of Eq. (|I.1[) are only miserably constrained. All available neutrino data allow one to "measure" 
M £ [0, 10"^] eV U M e [IQ-^, 10^-'] GeV, while the largest A^i e [2 x IQ-^^, 10]. It is clear that we have only begun 
to properly probe the origin of neutrino masses. 

In Sec. [n] we discussed a very useful basis in which to describe pseudo-Dirac neutrinos. In the basis of choice, 
the relation between the mass-squared splittings of the different quasi-degenerate states is directly related to the 
"diagonal" elements of Mij, while the "off-diagonal" ones are only responsible, at leading order, for affecting the 
neutrino mixing matrix. Bounds on these off-diagonal elements were not discussed in any detail, but we argued that 
these are much weaker than the bounds on the diagonal elements. It is worthwhile to investigate how such effects 
might be probed experimentally. 
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APPENDIX: MODEL INDEPENDENT FIT TO SOLAR NEUTRINO DATA, INCLUDING STERILE 

NEUTRINOS 

In this appendix, we discuss a model independent, simplified fit to the solar neutrino data, along the lines of analyses 
performed in the past by Barger, Marfatia and Whisnant [i^, H^]. In a nutshell, we "split" the oscillation probabilities 
of solar neutrinos into three disjoined regions, qualitatively constrained by different experiments: 

• P^: this is the oscillation probability of high energy solar neutrinos {E > 4 MeV), averaged over the entire 
energy region. In practice, is dominated by the ^B neutrinos, and is, by far, best constrained by the real 
time experiments SuperKamiokande and SNO. 

• P^: this is the oscillation probability of medium energy solar neutrinos (0.5 < E < 4 MeV). In practice, 
constraints on this region are dominated by the monochromatic ""Be neutrinos [E = 0.863 MeV). P^^ is best 
constrained by the Borexino experiment. 

• P^: this is the oscillation probability of the low energy solar neutrinos [E < 0.5 MeV). These consist almost 
exclusively of pp neutrinos and are only constrained by the Gallium experiments. 

Our analysis is very similar to that of ji^ HJ, to which we refer for details, with two exceptions. One is that the 
Borexino data was not available when [47[ was published and we discuss their impact here. The other is that we 
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FIG. 7: Allowed region of the P^^ x P^^ parameter space after combining data from SuperKamiokande and SNO experiments 
at the 1, 2, 3, and 4 a levels. Confidence levels are defined as regions of constant Ax^ = 2.3, 6.18, 11.83, 17.95 



wish to consider that the eleetron neutrinos can oscillate into sterile neutrinos, and discuss how constrained is this 
possibihty. 

Figure [7] depicts the allowed region of the P^^ x P^^ parameter space at different confidence levels after combining 
data from SuperKamiokande [111 and SNO [23|. From the SNO experiment we learn ^'^'^ = (1.67±0.09) x 10*' cm'^s"^ 
(j)^c: ^ (5 54 _^ 0.48) x 10^ cm-^s"!, while from SuperKamiokande we learn cj)'^^ = (2.35 ± 0.08) x 10^^ cm-^s"!. 
For the Standard Solar Model expectation we assume cj) = (5.69 ± 0.91) x 10^ cm~^s~^. We picked the theoretical 
prediction that leads to the most conservative bound on P^^ and assumed P^^ + P^^ + P^^ = 1 . After marginalizing 
over values, we obtain G [0, 0.37] at the three sigma level. 

Figure [5] depicts the allowed region of the P^j x Pj^/ parameter space at different confidence levels after combining 
data from Borexino [3l| and the Homestake experiments [s^. While computing the expected flux at Homestake, we 
use the value of P^^ obtained in the analysis depicted in Fig. [7] in order to estimate the contribution of neutrinos to 
the neutrino fiux measured at Homestake. For more details regarding how to treat the Chlorine data, we refer readers 
to [1^, dj. For the Borexino experiment, we learn that 49 ± 5 counts/day/100 tons have been observed, and that 
75 ±4 counts/day/100 tons were expected according to the Standard Solar Model. The constraint Pel + Pes + Pea = 1 
"cuts off" the part of the parameter space where Pj)f + Pes > 1. Note that our confidence level curves are defined 
as constant Ax^ = X^iPee ' Pes) ~ Xmin contours. One should hence be careful when translating what we refer to as 
an N sigma {N = 1, 2, 3, . . .) bound into a probability that a certain value is allowed. Marginalizing over Pjjf , we get 
Pes S [Oj 0.57] at the three sigma confidence level. 

It is interesting to note that, similar to the neutral current measurement at SNO and the elastic scattering result 
from SuperKamiokande, the Borexino result is sensitive to both electron-neutrinos and other active {Vfi^r) neutrino 
species, but Borexino data alone cannot rule out Pea = 0. In principle, the combination of Borexino and Homestake 
data is sensitive to whether Pj)f 7^ 0. Alas, this is not the case in practice, as one can readily see in Fig. [S] 

As noted above, all information on low-energy neutrinos comes from the Gallium experiments, and these are only 
sensitive to Pee (i.e. the Gallium experiments are only sensitive to electron neutrinos from the sun). On the other 
hand, since the Gallium experiments are sensitive to neutrinos of all energies above the experimental threshold, the 
Gallium result for P^^ depends on P^ and P^. P^ is well constrained by the *B experiments, so we concentrate on 
the interplay between P^^, P^^ and Pes (note that there is no experimental information on P^^ or Pea)- Fig. [9] depicts 
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FIG. 8: Allowed region of the P^^ x P^^ parameter space after combining data from Borexino and the Homestake experiments 
at the 1, 2, 3, and 4 a levels. Confidence levels are defined as regions of constant Ax^ = 2.3, 6.18, 11.83, 17.95. 



the allowed region of the x Pj|f parameter space after combining data from the Borexino and Gallium experiments 
(the contribution from Homestake data is also included. Its effect is very small). In order to illustrate the effect of 
sterile neutrinos, the figure depicts the one sigma allowed region obtained once we marginalize over all allowed values 
of PgY (boundary of the black region in Fig. [9|) and the two, three and four sigma allowed regions which are obtained 
if PgY is set to zero in the analysis. 

Our result is consistent with that shown by the Borexino Collaboration in the Neutrino 2008 conference [i^. It 
also depicts the correlation between P^^ and P^l , courtesy of the Gallium experiments. Finally the presence of sterile 
neutrinos in the fit renders the allowed region in the P^^ x Pj*/ plane larger, but does not qualitatively invalidate the 
results obtained assuming that there are no sterile neutrinos. 
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